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Abstract 

We study the correlators for interacting quantum field theory in the flat chart of de Sitter space 
at all orders in perturbation. The correlators are calculated in the in-in formalism which are often 
applied to the calculations in the cosmological perturbation. It is shown that these correlators are 
de Sitter invariant. They are compared with the correlators calculated based on the Euclidean field 
theory. We then find that these two correlators are identical. This correspondence has been already 
shown graph by graph but we give an alternative proof of it by direct calculation. 
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1 Introduction 

In recent years, there have been rapid progresses in the precise measurement of the observable quantities 
in cosmology, e.g., the non-Gaussianity of the fluctuations generated during inflation, which is expected 
to be a powerful tool as a probe of the early universe. Along with the development of these precise 
measurements, the need arises for the accurate theoretical predictions of the corresponding quantities. 

When computing the non-Gaussianity, one needs to discuss interacting quantum field theory on in- 
flationary background, in which one does not generally know how to define the interacting vacuum. One 
often uses the it prescription in cosmology to calculate the correlators perturbatively. (See, for example, 
Ref. pp.) In the Minkowski, this prescription is known to perturbatively give the Poincare invariant 
correlators for interacting theory, defining interacting vacuum as the lowest energy eigen state. Indeed, 
this prescription also enables us to calculate the non-Gaussianity or higher correlations in the inflationary 
era, but the physical meaning of it is not as clear as in the Minkowski case. Our main interest in this 
paper is in the meaning of the it prescription for interacting field theory in dc Sitter space. 

The free scalar quantum field theory in dc Sitter space is well understood [2H5] , while interacting one 
is a hot subject with a lot of debate We focus in the present paper on the problem whether the 

it prescription for interacting theory breaks de Sitter invariance. 

Since de Sitter space is maximally symmetric and possesses 15*0(4, 1), de Sitter symmetry, it is strongly 
expected to exist a de Sitter invariant vacuum even for interacting theory. In fact, a de Sitter invariant 
vacuum for interacting theory is defined by constructing arbitrary correlators perturbatively at all orders 
using the Euclidean method [TU]. While the vacuum state thus constructed are manifestly de Sitter 
invariant, it is not obvious whether the ones defined by the it prescription in the flat chart are de Sitter 
invariant. Notice, for example, that in the latter the integration region for the vertices in calculating 
correlators are restricted to the future of the cosmological horizon, which is not de Sitter invariant. 

Actually, this problem has been already resolved affirmatively in Ref. [TT] for interacting massive 
scalar field. Namely, the it prescription does not break de Sitter invariance for interacting massive scalar 
field. Furthermore, the vacuum defined by the it prescription has been shown to be equivalent to the 
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Euclidean vacuum. The main ideas in Ref. [TT] are as follows. They start from correlators defined on an 
Euclidean sphere and take, on the Euclidean sphere, coordinates such that when we Wick rotate the time 
coordinate continues to the static chart of the Lorentzian de Sitter space. Then, after the deformation of 
the integral path of the Euclidean time, fall-off of the propagator in the large separation limit leads to the 
identity of the two correlators at least on the static chart. From the analyticity of the in-in correlators 
for their time coordinates, and the uniqueness of the analytic continuation, it is shown that the in-in 
correlators in the fiat chart arc identical to the analytic continuations of those on an Euclidean shperc. 

Then, it is natural to ask what becomes of in massless field theory What happens for graviton in 
de Sitter space has especially been a topic of much discussion. (See, e.g. [T51l21fB5] .) Our final goal is 
to extend the correspondence between the two vacua to those interacting massless field theory. It is also 
worth considering derivatively interacting massless scalar field, which can be a step toward graviton. 

It seems difficult to extend the discussion of massive field theory above to massless field theory where 
the propagator does not fall off in general, since the proof of the correspondence between the two vacua 
relies on this decay property of the propagator at a large separation as explained above. In order to 
attack those theories, we take another approach. That is, we directly calculate the correlators with the 
ie prescription. We derive, along this way, the analytic Mellin-Barnes formulae for the correlators of 
quantum fields in the flat chart. The resulting correlators are shown to be completely the same as the 
analytic continuations of the ones considered in the Euclidean field theory in Ref. [TU] • Thus we find that 
the ie prescription in de Sitter space gives the vacuum state corresponding to the Euclidean field theory. 
Although we consider only massive theory in the present paper, we believe that our proof has potential 
to be extended to wider range of theories which include interacting massless theory such as derivatively 
interacting one, since it does not employ the decay property of the propagator. 

This paper is organized as follows. In Sec. [51 we briefly review how to describe de Sitter space, 
especially the flat chart, and massive free scalar quantum field theory on it. Pauli-Villars regularization 
scheme is also introduced. Then we proceed to the interacting theory, in Sec. [3J 0] and [SJ We consider, 
in Sec. [3] and HI a tree graph which contributes to an iV-pt correlator with single vertex. Then in Sec.[5j 
we extend the discussion to arbitrary graphs. We give a brief summary in Sec. [SI 

2 Preliminaries 

In this section, we briefly review free scalar quantum field theory on de Sitter space, especially in the fiat 
chart. We also introduce Pauli-Villars regularization scheme for later use. 

2.1 de Sitter space 

We consider D-dimcnsional de Sitter space dS D with, for simplicity, unit radius. This is a hypcrboloid 
embedded in (D + l)-dimensional Minkowski space with metric r\ a b = (— , +, • • • , +)• The embedding is 
specified by 

r, ab X a X h = 1 . (2.1) 

It is convenient to define the invariant distance between two points X and Y in de Sitter space by the 
Minkowski inner product of X and Y, which we denote as 

Z{X,Y):=r hb X a Y b , (2.2) 

as in Ref. [9]. For brevity, we often use alternative notation Zxy for Z(X,Y), Z\y for Z(Xi,Y) and so 
forth in the following. 
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The coordinates (77, x) in the flat chart are related to the embedding coordinates as 



x ° = \{«- 1 

2 V V 



2rj 2 V V 



2?y 77 



(2.3) 



where means the norm of (D — l)-vector x. The flat chart coordinates with —00 < 77 < and 
x £ R 15-1 span just a half of the whole spacetime region. In fact, the linear combination 



X° + X D 



1 

V 



is restricted to the positive side for negative 77. The metric in the flat chart is expressed as 

ds 2 = \i-dij 2 + dx 2 ) . 



(2.4) 



(2.5) 



Expressed in the flat chart coordinates, the invariant distance between X and X' , Z(X,X') is given by 

(v-rf) 2 - \\x-x'\\ 2 



Z(X,X') = 1 + 



2777/' 



(2.6) 



where (77, x) and (77', x') are the flat chart coordinates corresponding to X and X', respectively. 
2.2 Free QFT on de Sitter 

We now consider a massive free scalar QFT on de Sitter space. We focus on the Green's function G(X, Y) 
given by 



r{-a)r{a + D-l) ( D 1 + Z XY \ 



(2.7) 



which corresponds to taking Bunch-Davies vacuum |44j or Euclidean vacuum |45j . a is related to the 
mass of the field m by 



D - I 



D- I 



2 V V 2 

Expressing the hypcrgeometric function in the Barnes representation, we have 

1 — Zxy \ " 



G(X, Y) 



■ 



with 



■r 



_D 

2 



Here 



r 



T +(T, I" -7T -<7 



ai, a 2 , ■ ■ ■ 



(2.8) 



(2.9) 



(2.10) 



stands for r(ai)T(a2) ■ ■ ■ / T (/3i) T (fa) ■ ■ ■ , and the symbol /(■■■) means the Barnes integral. The 
Barnes integral is an integral along a straight line, C, that traverses from —zoo to +ioo parallel to the 
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imaginary axis with the factor \/2m: 



(-)-££(■■■). (-» 

The integrand of the Barnes integral includes sequences of poles. For example, r(z) possesses a sequence 
of poles at z = 0, — 1, — 2, • • • . The integration path C is taken to avoid all the sequences of poles in the 
integrand. In the case of the above Green's function, C is taken to satisfy 

max{-Rea-D + l, Recr} < Rei/ < minjl- y, oj . (2.12) 

This region of the integration path is called "fundamental strip," and the poles such that are associated 
with Gamma functions like r(- •• — !/) (/"(• • • + v)) and hence such that line up on the right (left) hand 
side of this strip are called right (left) poles. (See Fig. [TJ) The symbol like J is used to represent the 
Barnes integral in this meaning in the following. 



2.3 Pauli-Villars Regularization 

Because we consider interacting theory in the present paper, we have to introduce some ultra-violet 
regularization scheme. We make use of the Pauli-Villars regularization. This scheme attaches some 
massive propagators, Gi(X,Y), defined in Eq. (|2.9[) with m replaced by the regulator mass Mj, to 
the original one, G(X,Y), so that we replace the original propagator in a graph with the regularized 
propagator 

G les (X, Y) := G(X, Y) + ^2 C i G l {X, Y) . (2.13) 

i 

The coefficients C, are chosen so that the regularized propagator G leg (X, Y) becomes finite in the coin- 
cidence limit Y — > X, which leads to the conditions 

J2Ci = -l, £QM? = 0, £<7 4 M? = 0, ••• . (2.14) 

i i i 

This regularization scheme affects the pole structure of ipiy) in (|2.9p . eliminating the first several right 
poles of t\)(v) which are responsible for the behaviour of the Green's function in the coincidence limit [TU] . 
The regularized Green's function is written as 

G^(X,Y) = jf 0—^ ) r{-vW°*(y) , (2.15) 

where we assume that ■i/> reg (^) is regularized to be analytic in the region 

Rc(T<Rci/<p, (2.16) 

with p a sufficiently large positive constant. (See Fig. [TJ) In the following sections, we drop, for simplicity, 
the symbols such as reg on G and ip. 



3 Interacting QFT: Single Vertex 

We now move on to the interacting theory. The interacting QFT in the flat chart of the Lorcntzian de 
Sitter space is discussed in the present and the succeeding sections. When we express the correlators 
in the wave number representation, we employ the it prescription to calculate the correlators for the 
interacting vacuum. This prescription regularizes the oscillatory behaviour of the Green's functions at 
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Figure 1: The left figure shows the pole structure for ip{v) which is not regularized. There are two 
series of left poles from v = a and v = —a — D — 1, and right poles from v = 1 — D/2. The right one 
shows the pole structure for ip rcs (i/) which is Pauli-Villars regularized. The shaded region represents the 
fundamental strip in each figure. 

infinity in time and makes the vertex integral converge. Although what we discuss in the present paper 
is the position space representation of the correlators, we also employ the it prescription to specify the 
interacting vacuum. 

In this section, we discuss perturbative calculations of a single vertex tree graph for the correlators. 
Then, we identify the problems to be solved to accomplish this calculation, which are solved in Sec. |U 
In Sec. [5] the results for single vertex tree graphs are extended to arbitrary graphs. 

3.1 Definition of the In-in Path 

Let us consider iV-pt Green's function. The contribution to iV-pt correlator at the lowest order in 
perturbation theory is given by 



In the in-in formalism with the ie prescription, the integration region Q for the vertex integral is specified 
as follows. 

We first introduce an 77-integration path P on the 77-plane, independently of the spatial coordinates 
y, defined as a curve which starts from — ooe~ le and ends at — ooe" as shown in Fig. [2] All the external 
points are also supposed to be placed along this path. In case of the wave number representation, this 
construction completes the definition of the in-in path on the 77-plane. If we take the 77-path along P, the 
integral converges with the integrand vanishing fast enough in the past. 

For the purpose of the present paper, it is more convenient to use the position space representation 
to compute the correlators. The vertex integrals involve the spatial integration, too. As a starting point, 
we set the region of the vertex integral fi to P x R- 0-1 . If we first carry out the spatial integration 
before temporal one, the integral would diverge because we then pick up the contributions from distant 
spacelike region. On the other hand, if we integrate first for the time variable and then for the spatial 
ones, the integral is convergent as we see in Sec. |3J This means that the integral over P x R- 0-1 is not 
well-defined as a multiple integral. 

To make the integral to be well-defined as a multiple integral, we modify the integral region by 
deforming the path of the 77- integral P There are branching points on the 77-planc, which correspond 
to the intersections with the light cones emanating from the external points. On the 77-plane for fixed y, 
G(Xi,Y) has the same structure of Riemann surface as that of (1 — Z i Y)" i , where v.i is some complex 




(3.1) 
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Figure 2: This figure shows the ?7-path P which is later deformed to P y . The dots represent the time 
coordinate rji of the external points and the crosses are the branching points corresponding to the light 
cones emanating from the external points. The dashed lines are the branch cuts. 



number, and 



where 



i - = (-v + Vi,+)(v -m,-) 

2 4(-T7)(-7jO 



rj h ± := ?7j ± ||£Ci - y|| , (i = 1, • • ■ , N - 1, N) . 
Namely, the integrand has the same structure of Riemann surface as that of 



N 



(_ r? )-P+E^)J| ( _ 7? + J7 . +) ^ (r? _ 7? ._ )i 



(3.2) 



(3.3) 



(3.4) 



The time integration is unchanged even if we deform the integration contour as long as it does not 
cross singularities of the integrand. Thus, we deform the contour P to P y such that the maximum value 
of the real part of rj on P y is equal to max^jRe r)i t —} + b where 6 is a small real positive constant. (See 
Fig- El) This deformation on the ?7-planc is significant when the spatial coordinates of the vertex is 
largely separated from those of relevant external points. To the contrary, when \\xi — y\\ is small for i 
that realizes the maximum among Re the modified contour P y is almost identical to the original one 
P. 

Using this P y , we define the integration region 



(3.5) 



in C x R B_1 . The result of the integral is the same as that is obtained by integrating first for time and 
then for space for the original integration region, but we emphasize that the integral over H is now a 
multiple integral. 



3.2 Problems to Be Solved in the Calculations 

Let us return to Eq. (|3.1|) . Inserting Eq. (|2.9[) into Eq. (|3.1|) . we have 



Vjv(X[, • • ■ , Xjy) 







f dV Y f • 


■J " 


in Jv! 


J V N 



N 



' N 

n 

,i=l 



1 - Zc 



(3.6) 
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Figure 3: This figure represents the deformed contour P y for fixed spatial coordinate y. The original 
path P is deformed as long as it does not cross the singularities. 



If we can exchange the order of the integrals, J n dVy and J ■ ■ ■ J v , we are led to calculate the following 
integral 



M{v\, ■ ■ ■ ,1>N-1, vn) 



dVy 



1-Zn 



J NY 



(3.7) 



The first problem is to calculate this integral. This quantity is shown to have an analytic Mcllin-Barncs 
representation in Sec. SI and hence if this exchange of the order of integration is allowed, V/v can be 
represented in an analytic Mellin-Barnes form. It is not trivial whether this exchange of the order of the 
integration is allowed or not. This is the second problems. The same problem arises also for arbitrary 
graphs as for the tree level graphs. We will extend our discussion to arbitrary graphs in Sec. [5] 



4 Computation of the Master Integral 

The goal of this section is to compute the master integral: 

ka, ^ (1-Z 1Y \ V1 fl-Z n: 
M(vi,--- ,v n , vn) = J dVy ( ^ ) "' 

where we have introduced n := N — 1 for convenience, 

d V d D -\j 

is the invariant volume, and Q is defined in Eq. (|3.5[) . 



1 - Z 



NY 



4.1 Generating Function for the Master Integral 

In order to evaluate the above expression (|4.ip , we introduce the following generating function 



A(cei,- ■■ ,a„) := / dV Y V]< 
J n \i=i 



1 - Z, 



(4.1) 



(4.2) 



(4.3) 



following Rcf. |10j . in which it was used to evaluate the master integral on an Euclidean sphere. Here 

Re A < , ax , • • • , a n > , a N := 1 , (4.4) 
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Figure 4: A figure representing the in-in path P y for the external points which lie on the real Lorentzian 
section and are mutually spacelike separated. The crosses represent the branching points and the dotted 
lines the branch cuts. 



are assumed. 

In this subsection we establish the relation between the generating function and the master integral. 
Formally, in the same way as in the Euclidean case discussed in Ref. [TU], the generating function (|4.3[) 
seems to be related to the master integral p.7[) also in the present case as follows: 
[Step 1.] We first apply Eq. (|A.1|) to the integrand of (|4.3[) to obtain 



A(a lt ■■■ , a n ) = dV Y jf ( ai r ■■ j K) u " 



where 



( - — Zix_^ fl — ) ( 4 5 ) 



un ■= A - ^2 u i ■ (4-6) 

i=i 

[Step 2.] Next, we exchange the order of the integration, J n dVy and f u ••• J u , to have 

xT [-ui, • • • , -u n , -un] M (ui, • • • ,u n , u N ) ■ (4.7) 

Thus, the Mellin transform of A gives A4. 

However, we have to prove that [Step 1.] and [Step 2.] are indeed possible, which is the goal of 
this subsection. In particular, [Step 2.] requires that the integral over O is a multiple integral. The 
convergence of the integral is rather obvious when we consider the corresponding integral over a compact 
Euclidean sphere, while it is not in the present case where the integration region is non-compact. In this 
subsection, we assume, for a technical reason, that the time coordinates of all external points lie on the 
real Lorentzian section, i.e. r/, € R_, Hi € R 15-1 , and furthermore, that any pairs of them are mutually 
spacelike separated. 

Since the definition of the in-in path described in Sec. 13.11 requires the external points to lie along 
the in-in path and therefore their time coordinates are complex in general, we need some explanations of 
the in-in path for this configuration. The path is defined on the 77-plane by taking the limit Imrji — > 
in P y introduced in Sec. 13.11 It seems that the path in this limit must, at least partly, lie on the 77-real 
axis. However, since the external points are mutually spacelike, the branch cuts, lying on the 77-real axis, 
do not cover the whole 77-real axis. Therefore, the limit can be taken without the pass P y crossing the 
branch cuts, and hence the in-in path in this limit is simply a contour going from —00 e - " to —00 e le as 
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Vi.- 



m. 



ViA 



Figure 5: The dot represents the time coordinate of a vertex on P y . The remainder of the summation 
of the arguments of two vectors relevant to the subscript i and that relevant to the subscript j gives 
| arg(l - Z iY ) - arg(l - Z jY ) | as in (|4T0)) . 



shown in Fig. [4] 

Proof of [Step 1.]: Note that the following inequalities hold for arbitrary Y € fi: 



arg(l - Z lY ) - arg(l - Z jY ) | < tt , 



1,2,- 



In fact, arg(l — Z{ Y ) is given by 

arg(l - Z iY ) = axg(-r) + + arg(?7 - % _) - axg(-ri) - arg(-%) 



(4.8) 



(4.9) 



and then, noticing that arg(— rji) = since all the external points are on the real Lorentzian section, we 
have 



arg(l - Z iY ) - arg(l - Z jY ) \ 

= I axg(-r? + + arg(ry - ^ _) - arg(-?y + - arg(r/ - %,_) | . 



(4.10) 



This quantity is less than ir for any (rj, y) € f2. (See Fig.0) The inequality (|4.8p is the sufficient condition 
that the formula (|A.1[) can be applied to the integrand of Eq. (|4.3[) . For the later purpose, we modify the 
integration path P y as such that satisfies 



arg(l - Z lY ) - arg(l - Z jY )\ < tt - S 



(4.11) 



for any i and j with a small positive number 5. This can be achieved easily. Because | arg(l — Z.;y) — 
arg(l — Zj Y )\ is close to tt only in the small region surrounding the interval {r)i.-, ?7j.-) or (^i,+ , %.+), the 
path can be chosen to avoid this region. 

Proof of [Step 2.]: We denote the integration paths for u\,- ■ ■ ,u n as C\, • • • , C„, respectively, and 
define C := C\ x • ■ ■ x C n , The sufficient condition to allow to exchange the order of the integration, 
j n dV Y and jcX\^ = idui/2TTi, is that the integral is absolutely convergent (Fubini's theorem). In the 
present case, we should examine the following integral 



\r(-x)\J c 



n 

fc=i 



du k 



2ni 



(ai) Ul •••(«„)""! \r[- Ul , ■■■ , -u n ,-u N }\ 



\dV Y \ 



i-z 1Y \ U1 



\-z NY \ UN 



(4.12) 
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where 



\dV Y \ 



\dr)\d D -\) 



n 



D 



(4.13) 



If this integral is finite, then we can justify the exchange of the order of integrals in [Step 2.]. 

To show this, we focus on the integrand of the Vt integral in the large brackets in Eq. (|4.12p for fixed 



1 - Z 1Y 



1 — Zmy \ ' 
2 ) 



Notice that 



(l-Z iY ) u *\ = \l-Z lY \ Rcu > cxp [- arg(l - Z iY )Im Ui 



(4.14) 



(4.15) 



Along the integration path of Ui parallel to the imaginary axis, Imitj varies while Rcji, is fixed. Taking 
into account that um includes Uj as given in Eq. (|4.6p . the part depending on Imu,; in Eq. (|4.14p is 
factored out as 

cxp[{arg(l-Z J vy)-arg(l-Z l y)}Imu ? ] . (4.16) 



Since |arg(l — Z^ Y ) — arg(l — Zfy)\ is bounded as shown in Eq. (|4.11jl . this factor is bounded from above 
by exp [(it — S)\Imui\] . Therefore, noticing that on is real positive number, wc find that 



[Eq. (I4T21) ] < 



l^(-A)| 



ai) RoMl 
du k 



(a n ) Rcu -\ I \dV Y \ 
In 



1 - Z 



1Y 



Rc ui 



1 - z 



NY 



Rc ujv 



n 

fe=i 



2ni 



\r[-ui, 



-U n , ~u N \\ e 



(7r-5)(|ImMiH h|Imu„|) 



(4.17) 



Since \r(x + iy)\ ~ (2n) 1 / 2 e n \ y \/ 2 \y\ x 1 / 2 (|y| — > +oo), u k integrals in the second line in the last expres- 
sion are convergent. Therefore, our remaining task is to show that the volume integral 



/ \dVy\ 



(i^) 


Re u\ 




Re u n 











(4.18) 



is also finite. 

For this purpose, we first introduce a representative point Xq with coordinates in the fiat chart defined 

by 



N 



(Vo,x Q ) := ^2pi(m,Xi) , fa > 0, = 1) , 



(4.19) 



and a domain Dq far from Xq in terms of the invariant distance by 

D :={Y\\Z 0Y \> Z Q }nn . (4.20) 
Note that if we take Zq to be sufficiently large, we see that 

\Z lY \ > const, x \Z 0Y \ , (Ye D , t = 1,2,--- ,N) . (4.21) 
We divide the region f2 into (tt\Do) and Dq, and evaluate each contribution to (|4.18[) separately 

(i) Integral over fl\Do: 
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Wc further divide £1\Dq into K denned by 



K := {(v,y) e Cl\D | ||aj - y\\ > R} , (4.22) 

and its complement (Q\Dq)\K . R is set large enough for K not to include any external points. (See 
Fig.©) 



(i-a) Integral over {V£\Dq)\K: 

The region (D,\Dq)\K is compact but it contains the coincidence points (rj, y) — (r)i, Xi) at which the 
integrand of (|4.18[) diverges. Since y « Xj around them, the path P y is identical to the original one P, 
and hence r\ = t]i + iSrj with real Srj. Then, we have 

|1 - Z iY \^ « ((<5 ?? ) 2 + ||asi - y|| 2 ) RGUl (4.23) 

around a point (r/j,x,), which shows that (|4.18p is finite as long as we choose the integration path of itj 
to satisfy 

Re Ul >-D/2, (i = 1,2,- •• ,n,JV) , (4.24) 

which does not conflict with [Step 1 .] . Recall that the fundamental strip of Eq. (|4. 5[) contains the paths 
with Re Mi for all i being infinitesimally small negative constants. 



(i-b) Integral over K: 

We first see that, for Y G K, |1 — Zjy| is bounded both from below and from above by positive 
constants. Recall that the 77-path P y is defined by deforming P not to touch r]i t - except for the case with 
y m Xi, which occurs in (Q\Dq)\K . Therefore, |1 — Ziy\ does not vanish, bounded from below by some 
constant c_( > 0). It is also easy to show that |1 — Z^y \ is bounded from above by some constant c+. If 
jl — Zfy\ is sufficiently large, Zqy will be larger than Zq. Then, by the definition of K, Y is not included 
in K. Thus, we conclude that for some positive constants c±, 



c_ < \l-Z tY \ <c+ , (YeK). 
Furthermore, one can claim that the volume of the region K is finite, i.e., 



/ \dV Y \ < +00 

J K 



(4.25) 



(4.26) 



In showing this, the non-trivial point is that the region K extends to infinitely large ||y||. However, the 
region of the 77-integral is confined to the interval 



b' < Re i] < r?i( y ),_ + b 



(4.27) 



where b is the same constant used in defining the path P y and i(y) is the label of the external point such 
that > Tjic,- for all k 7^ i{y)- Here the point is that one can choose a large positive constant b' to 

be independent of y. In fact, the invariant distance between Xq := (t}o,Xq), and the point corresponding 
to the above lower bound Ybdry := (wiu,\ - — b\y) = (rji(„\ — Wxi 



|l-Z(X ,F bdry )| 



y (%(»),- - b ',y) = (r}i( y ) - \\Xi( y ) — 2/|| — b',y) is evaluated as 
(Vo - Vi(y) + \\xi(y) ~y\\ + b') 2 - \\x i[y) - y\\ 2 



> 



m\ 



2 Vo(Vi{y) ~ \ \ x i{y) -V\\- b') 
(\\ x i(y) ~y\\ + b ' -> +00) 



In the last inequality we assumed ||a;i( w ) — y\ \ + b' 



(4.28) 

00 but this should be a good approximation in 
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Figure 6: This figure is a schematic of how we divide the integration region Q. There are -Do, K 
and (Q\Do)\K. The dots except for X represent the external points. The dashed lines represent 
schematically the "past light cone of Xq." 



the region K . Therefore, if b'/( \— r)o\) is taken sufficiently large compared with Zq, the above range of r\ 
covers the whole region of K. Thus, the volume J„ \dV Y \ is bounded by 



\dVy\ <ci 



K 



d D -h 



||aj -y||>fl 



\dn\_ 

nny),--v \^\ D 



< C 2 



<c -v\\>R W x i(y) y 



\D 



< + OO 



(4.29) 



where c\ and ci are some appropriately chosen constants of 0(1). In the second inequality we used 
\r]\ > \f)i(y),- + b\ ~ ||aJi(y) — y\\- Therefore, the integral over K is proven to be finite. 



(ii) Integral over Do: 

We next proceed to the integral over Dq. Using Eq. (|4.21[) . one can easily bound the volume integral 
of our current concern from above as 



/ \dVy\ 




Re ui 




J D 









Rc uj\ 



< c 3 x T 



\dV Y \\Z 0Y \ RcX 



(4.30) 



where C3 is a constant of 0(1), and we have used the relation Ui = A. 

In order to show this integral is finite, wc use Zqy as a time coordinate instead of rj, which leads the 
integration measure to transform as 



d V d D -\j= I 1- 



Va z 0Y 



r) d{Z 0Y )d D \j 



v/l^o-ylP + ^^oV- 1 ), 



(4.31) 



We substitute this into the right hand side of (|4.30[) . Approximating Zq Y — 1 ~ Zq Y and introducing 
x := y 1 1 {—tiqZqy), we find that the integral is finite as 



\dV Y \ \Z t 



0Y\ 



I Re A 



< C 4 



dZ 



Z 



gl—Re X 



dx ■ 



/y/l + x 2 \{l + ^l + x 2 ) 



D-l 



< +00 



(4.32) 



where again C4 is a constant of 0(1). 
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(iii) Summary: 

We have shown in this subsection that the integral (|4.12[) is indeed finite when the external points 
X\ , • • • , X n and Xn lie on the real Lorentzian section and are mutually in spacelike separation, as long 
as the integration contours for u\ , ■ ■ ■ ,u n satisfy the additional conditions (|4.4p and (|4.24l) : 

iV 

Re A = Re V] Ui < , 

i 

Reui>-y, ••• , Reu n > -—, Reu N > -— . (4.33) 

Then, the order of two integrals J n dVy and j c Y\^ =1 dui/2iri in Eq. (|4.5p are exchangeable, which 
implies that the master integral M. is given by the repeated Mellin transform of A. Furthermore, under 
these conditions A{a\, ■ • • , a n ) is finite and thus from Eq. (|4.7[) the master integral Ai(u\, ■ ■ ■ , «„, un) 
is also finite. That is, the master integral Ai(ui,--- ,u 71 ,un) is finite when the external points are in 
the real Lorentzian section and are mutually in spacelikc separation, with the conditions (|4.33l) satisfied. 
The analytic expression for Ai is given in the succeeding subsection, where the conditions on the external 
points are relaxed. 



4.2 Calculation of the Generating Function 

We now proceed to compute A and hence M. , to show its equivalence to the analytic continuation of the 
Euclidean correlators. Again in this subsection we first assume that all the external points Xi lie on the 
real Lorentzian section and that they are mutually in spacelike separation. After that, we show that the 
time coordinates of the external points r\i in the obtained expression for M. can be analytically continued 
to any point on the in-in path. 

The expression for A given in Eq. (|4.3|) can be transformed into 

**■■■■■*•> - / R „_/ D v,„T^(f>-H ■ <4 - 34 » 

where V ■ Y is an inner product of V and Y with respect to {D + l)-dimensional Minkowski metric and 

N 

V = Y,^X t . (4.35) 



i=l 



Notice that 



N N 

V° + V D = V ^ > , V = Y ^-xi . (4.36) 

' -iji j^i -m 



1=1 



Setting 



V ■ Y can be expressed as 



T/0 , yD N 

V, R:=(V + V D )x-V = J2—(y-^), (4-37) 

1 i=i ^ 



y. Y = . T + ^f±l, (4 .38) 
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where 



JY 



TV 



V -V := Vab V a V b = X) a * + 2 E° 4a J'( Z «>' - ^ 



(4.39) 



Thus, we obtain 



^4(ai, ■ • • , a„) 



.-A 



■„ (-»)' 



(tr- VY ) 



27ri x 2- I3 + 1 - A 



2 



dr(-r) 



— D 



E 

i=l 



a, + t 



(It 



pi 2iri 



-A-t)- u -\-t + t + )\t-t-)< 



(4.40) 



where P' y is the scale-transformed path of P y by a factor of (V + V )/2, and 

JV JV 

r± := -{-F± V-R 2 + J 2 } , F:=^a 4 , J 2 := 2 ]T ^(1 - Zy) 

i—1 i<Cj 



(4.41) 



As is mentioned at the beginning of this subsection, we have assumed that the external points are all 
mutually in spacelike separation, so that 



J 2 > , t± G R . 



Changing the integration variable further to £ := r — r_ , we obtain 



*4(ai, ■ • • , a„) = — 27ri x 2 



-D+l-A 



\\<-\ 



2m 



where 



.4 



i + \/i2 2 + J 2 } , B := r+ - r_ = v^ 2 + J 2 , 



(4.42) 



(4.43) 



(4.44) 



are both positive, and C is the integration contour shown in Fig. [7j which corresponds to P' y but the 
direction is reversed. (See Fig. [4] as a reference.) Notice that the path C does not have to respect the ie 
prescription here, because this integral is convergent without relying on the ie regulator. 
In order to compute this integral, we use the formula 



r 



-a + ii, -(a + (3 + 7 + 1) + fx, -/i, a + 7 + 1 — Li 
-a, -(3, -7, 7 + 1 



which is valid for Re (a + f3 + 7 + 1) < 0. 

It is not difficult to verify (|4.45[) . We denote this integral as I (a, f3, 7). When Re 7 > —1, the path 
C can be contracted to the forward and backward paths along the negative real axis. Noticing that only 
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Figure 7: The path C in the integral (|P5|) . 



the argument of £ 7 changes between these two paths, one can transform I(a, /3, 7) as 



I(a, /3, 7) 



2ni 





+00 

1 



{-dx)e- l ^x"t{x + A) a (x + Bf + 



2ni 



(-dx)e l ^x~< 



r[- 7 ,7 + l] Jo 



dx (x + A) a (x + B)^x' y 



Next, we expand (a; + 5)^, using Eq. (|A.2|1 . as 



{x + Bf 







'ft 


-0 



c^B 13 -* 1 , (Re/3 < Rc/x < 0) . 



Substituting this into Eq. (|4.46[) , we carry out a; integral first to obtain 



/(a, /3, 7) = / r 



-fS + fj,, -fi, 7 + 1 + fi, -a - 7 - 1 - fi 
-7, 7 + I, -{3, -a 



(4.46) 



(4.47) 



(4.48) 



Of course, the convergence of x integral imposes a condition Rc (a + 7 + 1 + fi) < 0. This is in fact 
satisfied because we can set Re/x arbitrarily close to Re/3 as long as Re/3 < Re/x(< 0) is maintained. If 
we change the integration variable from fi—¥fi — a — 7 — 1, we obtain the expression (|4.45[) . 

Finally, we remove the restriction Re 7 > — 1. In fact, the integrand is analytic for 7 and the £ 
integration is uniformly convergent for 7, as long as Re (a + j3 + 7 + 1) < 0. Therefore, the integral is 
analytic for 7, which enables us to remove the restriction RC7 > —1 by analytic continuation. 

Substituting Eq. (|Q5jl into Eq. (|Q5]l . we find 



-2m x 2-°- x+1 I d u - l R I r 

ft 



{R z + J 2 ) = 



A(ai, • • • , a n ) 

x2~" \F+ VR 2 + J 2 
We next carry out R integration, using the formula 

d D ~ l R [R 2 + J 2 ] v ' 2 \F + VR 2 + J 2 

r 



D + A + fi, D - A - 1 + /x, -/x, —D + 1 - fi 
D + X. -A. -A. A + l 



(4.49) 



R D -i 



-fi + K, -K, 



-I. + 

2 



s (4,50) 

which is valid when Re (fi + v + D — 1) < 0. The idea of the proof of the above formula is not so different 
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from that of the formula (|4.45|) . One applies Eq. (|A.2|) to F + \/ R 2 + J 2 to obtain 

F"- R (R 2 + J 2 ) K / 2 , (Re^<Re K <0). (4.51) 



F + v 7 /? 2 + J 2 


'-I r 


— /i + K, —K 






J K 







Substituting this into the left hand side of Eq. (|4.50[) , we obtain 

d D - l R [r 2 + j 2 ] v ' 2 If + Vr 2 + j 2 } " 





— jl + K, — K 


pH — k ju+n+D — 1 x 


1 K 







d~Z^(l + Z)— ; (4.52) 



where we have introduced a new integration variable S := (||i?||/J) 2 , and 



D-2 = 



2tt— 



(4.53) 



is the surface area of D — 2 dimensional unit sphere. The 3 integral is convergent if Re (v + k + D — 1) < 0, 
which can be satisfied since we can choose ReK arbitrarily close to Re/i as long as Rep, < Re«;(< 0) is 
maintained. Integration over E leads to (|4.5QD . 

Applying the formula (|4.50[) to the expression for the generating function (|4.49[) . and replacing the 
integration variables k and p. respectively, with 



k — /i + A 

w := and p := p, + D — 1 



(4.54) 



2w — A, — w 



D-l 
2 ' 



we obtain 

A( ai ,--- ,a n ) = (-i)2 2 - A 7T^ i / r 

J w 

x / 2~ p r [-A + p, A + 1 + p, -p, A + D - 1 - 2w - p] 
J P 

Finally, we perform the p integration in the above expression for A, using the formula 



2A+a-2 

2- p r[\ + 1 + p, -A + p, -p, A + a - 1 - p\ = — 

\/1T 



-A.A + 1,— ,A + - 



(4.55) 



(4.56) 



which can be proven as follows. If we close the p-path on the left hand side of (|4.56[) to the right, we 
have 

/ 2- p T[X + 1 + p, -A + p, -p, A + a - 1 - p] 
J P 

= r[-A, A + 1, A + a - l] 2 Fi ( -A, A+l;2-A-a;^ 



,1-A- 



a F[2A + a, o - 1, 1 - A - a] 2 F 1 2A + a, a - 1; A + a; - 



(4.57) 



Now applying the following formulae, known respectively as Bailey's summation theorem and the Gauss' 
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second summation theorem [46], 



2 F 1 ( a, 1 - a; 7; r ) = 2^^ 



7+g 7+(l — a) 
2 ' 2 



2a, 2/3; a + /3 + i; i ) 



we obtain after simple calculations (|4.56[) . Substituting a = D — 2w in (|4.56[l , we find 



A(a u --- ,a n ) = {-i)(^) D ' 2 / r 



2w — A, — w, A + — — w 
D + A, -A 



F 



X-2-u 



J 



(4.58) 
(4.59) 

(4.60) 



Recalling the definition of F and J, (|4.41j) . we expand F x ~ 2w J 2w to be integrals with respect to 
the power law indices of ojj's using Eq. (|A.1[) [10) . Since the master integral M. is given by the Mcllin 
transform of the generating function A, we finally obtain 



-i) 



(47T) 



D/2 



D 



j-j- / 1 



r(-%) 



[n^( 



r u+E^-E^ 



(4.61) 



where ^ represents AT (AT — l)/2-hold integration 



<i<j<N 



(/ dhij/2ir^J, and 



Zf< := ^ /i fcl + 51 



(4.62) 



fc=i+i 



In the above derivation of the equivalence between the expressions (|4.1[) and (|4.61|) we assumed that 
all external points are mutually in spacelike separation. However, we can easily extend the result (14.611) 
to the case of timelikc separation. First, notice that the integrand of Eq. (|4.1[) is analytic for the time 
coordinates of the external points r)i and also that this fi integral (|4.1[) continues to be well-defined and 
uniformly convergent even if rji are analytically continued to the region of timelike separation. On the 
other hand, the (/ly) integrals in (|4.61|) are convergent as long as | arg(l — Zij)\ < ir. This condition is 
satisfied when rji are placed on the original path P. Later, we need to replace Xj to Y and then Y is 
placed on the path P y . Even in this case it can be easily verified that the conditions | arg(l — )| < 7T arc 



satisfied. As a result, by the uniqueness of the analytic continuation, the master integral (|4.1[) is identical 
to the expression (|4.61[) even if the separations between some pairs of external points are timelike. As 
a remark already mentioned at the end of the preceding subsection, the parameters u\, ■ ■ ■ ,v n , must 
satisfy the conditions (|4.33|) in order for M.{y\, • • • , v n , vn) to be defined. 

Furthermore, without violating the convergence conditions, we can continue the external points in 
(|4.01l) to the Euclidean region where arg(l — Zij) = 0. Then, we find that the expression (|4.61[) is 
identical to the one obtained in the Euclidean Field Theory in Ref. [TU], except for the factor of — i due 
to convention. 



5 Interacting QFT: Arbitrary Graphs 

In the preceding section, we computed the master integral for a massive scalar field using the in-in for- 
malism in the Lorentzian dc Sitter space with the it prescription assuming the Euclidean vacuum at the 
level of non-interacting theory. We found that the resulting master integral is the analytic continuation 
of the one computed by the Euclidean path integral. Then, it might be expected that these two pertur- 
bative correlators are equivalent to all orders of perturbation. We will prove this equivalence along our 
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formulation in this section. Note that this equivalence is already shown to all orders, graph by graph in 
Ref. in a strictly different way from the present paper. 

5.1 Statement to be proven by induction 

It is known that the Euclidean path integral gives us a certain analytic form corresponding to any graph 
Vn(X±, ■ ■ ■ ,Xn), which contributes to the iV-pt correlator. The analytic expression for Vn is found in 
Ref. [10] in the form 

V N (hij), (5.1) 

where V^ihij) satisfies the following properties: 

1. The fundamental strip for each variable hij of V^{hij) contains the region 

Rehij G (a — V%j(h'), 0] , (5.2) 
where Vij is a linear combination of Re hki excluding Re hij with non- negative coefficients . 

2. When hij is in the region (|5.2p . Vj\r(/ijj) falls off, for fixed hm except for hij, as rapidly as 

V N (--- , hj =x + iy, ■■■)^e-^/ 2 \y\ x - 1 (\y\ > 1). (5.3) 

In this section we shall show by induction that any correlators calculated in the in-in formalism have the 
same analytic form as the above obtained in the Euclidean path integral. 

We start with some (N + K)-pt correlator Vn+k(X\, ■ ■ ■ ,Xm+k) which satisfies the properties 1 and 
2 above. The succeeding steps are as follows: 

(a) Set K external points, -Xjv+i, ■ ■ ■ ,Xn + k, in Vn+k to Y. 

(b) Add M — N propagators connected to Y and integrate over f2 with respect to Y, which gives a new 
M-point correlator with more loops. 

Any graphs can be obtained by this construction, except for the ones containing "one-link" loops, 
which are to be renormalized. It has been already shown in Ref. [10] that the intermediate (N + l)-pt 
function obtained in step (a) satisfies the properties 1 and 2. Therefore, what we have to consider is step 
(b). The resulting correlator, Vm(X\, • • • ,Xm), is given by 

Vm{X x ,~- ,X m ) = f dV Y V N+1 (X lr -- ,X N ,Y)G(X N+1 ,Y)---G(X M ,Y). (5.4) 

Integration region Q is specified in the same manner as in Sec. 13.11 but now with M external points, 
X\, • • • , Xm- We show below that the M-pt correlator given in Eq. (|5.4[) has the form of Eq. (|5.ip . 

5.2 Proof 

We here set the external points in Eq. (|5.4[) . X\,- ■ ■ ,Xm, to lie on the real Lorentzian section with 
mutually spacclike separation for technical reasons as in Sec. [4] Once we succeed in proving that 
Vm(Xi, ■ ■ ■ ,Xm) in Eq. (|5.4j) satisfies the properties 1 and 2, it is obvious that the time coordinates 
r\i (7 = 1, • • • , M) in Va/(Ai, • • • , Xm) can be analytically continued to the timelike separation or the 
Euclidean region for the same reason as we discussed for M. in the preceding section. 

1 In Ref. 1101 . Vij is set to be "a polynomial function of all Re h^i except for Re hij with non-negative coefficients," which 
does not matter here. 



Vn[Xi, ■ ■ • , Xn) 



JV 

n 

i<j 



1 - z. 



r{~h l3 ) 
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Representing the respective factors in (|5.4[) in the Mellin-Barnes form, i.e., (|2.9[) for G's and (|5.ip for 
Vjv+i, we obtain 



V M 



dV Y 

J(hij) J[vi] 
M 



n 



N , fy \ hi 

l-Zr Y 



n-hij) 



' N 

n 



i' J U'=iV+l 



r(-v r )ij)(vr) 



1 - 



VAr+i(/ijj,^) 



(5.5) 



where we have set the variables in the Barnes integral of G(Xp,Y)'s (I' = N + 1, • • • ,M) to uji and 
those for Vjv+i to hij (1 < i < j < N + 1), and we replaced /i,,jv+i (1 < « < JV) with i^. Here, we have 
denoted in short, 



(•••):= / •••/ (•••) 



(5.6) 



and so forth. The integrals for hij and Kj is a multiple integral and here we refer to the integration region 
for them as C. We rewrite Vm above, using C, as 



V M = 



ndhij -pi- dvi 
J ^i<j t=l 



1 - 



n 



r(-^) 



.7=1 



VN+i{hij,i>i) 



M 

n ^(^o • 

I'=N+1 

(5.7) 



Now the question is whether Q integration and C integration are exchangeable. In order to examine 
it, we take the absolute value of the integrand and repeatedly integrate it to see whether the integral is 
finite or not. Here we consider the following repeated integral 



II - II 



Wy\ 



2tti 

n 



n-hj) 



\V N+1 {h lJlVl )\ [JJ|V»( Vl')\ 



(5.8) 



where we have dropped the indices for Y[ and J^. As a default, the ranges of various indices are understood 
as 1 < i < N, N + 1 < V < M and 1 < I < M. 

In order to evaluate f2 integration of the above expression, we apply the discussion in Sec. 14.11 but 
slightly modify it. In Sec. 14.11 the essential point is the bound for | arg(l — Ziy) — arg(l — Zny)\ because, 
in Sec. 14.11 the exponents, itj, of (1 — Z;y)'s in the integrand are not independent since ^ Ui = A. 

In this subsection, however, they are mutually independent. Therefore, we should evaluate | arg(l — 
Zjy)\ itself. In fact, noting that the part dependent on the external points of the integrand of the H 
integration is expressed as 



J] |(1 - Z IY )\ Rcu ' exp [-5>rg(l - Z IY )lmui 



(5.9) 



we have to bound | arg(l — Ziy)\ for our purpose. 

Since the regions on the 77-plane where | arg(l — Zjy)\ = tt arc half lines going from —00 to r//._ or 
from r)i t + to +00 on the real axis, | arg(l — Ziy)\ is obviously less than 7r if Y e fl. Furthermore, since 
the path P y on the 77-plane is, by definition, tilted by e in the far past and deviates finitely from the 
region above, | arg(l — Ziy)\ is bounded by 7r — 5' with S' being some finite positive constant. Therefore, 
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we can factor out the Im vi dependent part in Eq. (|5.9[) to obtain the bound 



[Eq. (ED] < [] 1(1 - Z IY )\ KCUI exp [(tt - <50]T|W/I 



(5.10) 



Now the convergence of the volume integral follows in the exactly same manner as before. Thus, we are 
led to discuss the following integral 



dvi 



2ni 



n 



1 - Zi. 



[n^(^)i 



cxp 



n-hij) 

(7r-<5')^|Im^| |r[-i^,... ,-u N }\ 



x|Viv+i(/H7,fi)| 
Recall that ^(v) behaves as 

Mx + i V )\ -> e-^y^ivr 1 , (|y|»l) , 

and \r(x + iy)\ « ^Tr) 1 / 2 ^*"^/ 2 !?/!* -1 / 2 (|y| -> +oo). Furthermore, V N+1 behaves as 
\V N+1 (h 12 =x + iy, • • • )| -> e-'^'^lwr- 1 , (|y| » 1) , 



(5.11) 



(5.12) 



(5.13) 



from the property 2 of the assumption of induction, and the same is true for the other arguments, too. 
Therefore, the integral (|5.11j) is convergent, and hence the order of the integration over f2 and C in (|5.7[) 
are exchangeable: 



J J 



n 



'{ha) J[ui] 

X [ dVy 



r{~h t] ) 



M 

n 



l-Z IY 



(5.14) 



Now substituting the Mellin-Barnes form for the master integral, (|4.61[) . into Eq. ()5.14p . we arrive at 
the same Mellin-Barnes representation for Vm as that obtained in Sec. 4.2 of Ref. [10] , where Vm is shown 
to be represented in the form of Eq. (|5.1[) with Vm satisfying the properties 1 and 2. This completes the 
proof of the equivalence between the two types of correlators. 



6 Summary 

In this work, we considered massive interacting scalar field theory and demonstrated pcrturbativc calcu- 
lation for the correlators using in-in formalism in the flat chart of de Sitter space with the ie prescription. 
We found that the master integral defined in Eq. Q3.7P has completely the same Mellin-Barnes represen- 
tation as that obtained in Ref. [TU] based on the Euclidean field theory. We then derived the analytic 
Mellin-Barnes formulae for the correlators of quantum field on the flat chart. The resulting correlators 
are shown to be completely the same as the analytic continuations of the ones considered in the Euclidean 
field theory. Thus we find that the ie prescription in de Sitter space gives the Euclidean vacuum. 

Although the relation between these two vacua has been clarified in Ref. |llj . in order to extend 
this to massless field theory, we gave an alternative proof of their equivalence by direct calculation. In 
particular, graviton in de Sitter space has been a topic of much discussion. (See, e.g. [T3"ll2"lTf2"3] .) It is 
also worth considering derivatively interacting massless scalar field as a model of graviton. 

The proof in Ref. [11] of the equivalence between the two vacua relies on the decay of the propagator 
at a large separation. But the propagators in massless theory do not fall off in general. This could be 
an obstacle in extending the discussion to interacting massless field theory. Though we considered only 
massive theory in this work, we believe that our proof has potential to be extended to wider range of 
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theories which include derivatively interacting massless field theory, since our proof of the correspondence 
of the correlators is based on direct calculation without relying on this property. 
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A Formula 

Let A\,--- ,A n+ i be complex numbers satisfying |argAj — axg Aj \ < tt (Vi,j). Then, the following 
formula is true as a repeated integral and also as a multiple integral since the integral is easily shown to 
be independent of the order of the integration: 

{Ai + A 2 + ••• + A n+1 ) x 

= TF\)J U r [- A + E u « -«!."■ >«n] (A 1 ) Ul ---(A n ) u "(A n+1 ) x -^ u \(AA) 

Proof of (|A~Tj) : 

The basic formula is the following: 

(a + b) x = j^ J r[-X + fx, -^b x -^ , (|arg a -arg&| < tt) . (A.2) 

One applies this formula (|A.2[) with a = A n , b = A\ H — • + A„_i + A n+ i, and then again apply (| A.2|) to 
(Ai + - ■ • + A„_i+A n+ i) A ~ Al in the result of the previous step with a = A„_x, b = Ai + - ■ ■ + A„_ 2 + A n+1 . 
Repeating the same operation, one formally reaches (|A.1[) . The point is that the conditions 

| argAi - argA„ + i | < tt , 

| arg A 2 - arg(Ai + A n+1 ) \<it , 

| arg A n - arg(Ai + • • • + A„_i + A n+1 ) \ < tt , (A.3) 

are required to perform the above transformation. To allow the exchange of the order of the repeated 
integration without changing the result, we impose stronger conditions 

I arg A P{1) - arg A n+1 \ < tt , 

I arg Ap( 2 ) - arg(A P(1) + A n+1 ) | < ix , 

I axgAp (n) - arg(A P(1) H h Ap {n _ 1} + A n+1 ) | < tt , (A.4) 

for any permutation P. It is easily verified that, if we choose Aj's to satisfy 

| arg A, - arg A j \ = | sxg{Ai/Aj)\ < ir/2 , (A.5) 

for all pairs of i and j, the conditions (|A.4[) are all satisfied because in general \sag(J2i r i)\ < 7r/2 if 
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I arg(rj)| < ir/2 for all rj where rj's are understood as Ai/A^s. However, this restriction can be easily 
relaxed by analytic continuation with respect to Ai as long as the conditions 

\&rg Ai - arg^-| < tt , (A.6) 

are satisfied for all pairs since the right hand side of Eq. (|A.1[) continues to converge under these conditions. 
This completes the proof of (|A.1[) . 
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